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Abstract. Realistic ab-initio 3D, radiative-hydrodynamical 
convection simulations of the solar granulation have been ap- 
plied to Fe I and Fe II line formation. In contrast to classical 
analyses based on hydrostatic ID model atmospheres the pro- 
cedure contains no adjustable free parameters but the treatment 
of the numerical viscosity in the construction of the 3D, time- 
dependent, inhomogeneous model atmosphere and the elemen- 
tal abundance in the 3D spectral synthesis. However, the nu- 
merical viscosity is introduced purely for numerical stability 
purposes and is determined from standard hydrodynamical test 
cases with no adjustments allowed to improve the agreement 
with the observational constraints from the solar granulation. 

The non-thermal line broadening is mainly provided by 
the Doppler shifts arising from the convective flows in the so- 
lar photosphere and the solar oscillations. The almost perfect 
agreement between the predicted temporally and spatially av- 
eraged line profiles for weak Fe lines with the observed pro- 
files and the absence of trends in derived abundances with line 
strengths, seem to imply that the micro- and macroturbulence 
concepts are obsolete in these 3D analyses. Furthermore, the 
theoretical line asymmetries and shifts show a very satisfac- 
tory agreement with observations with an accuracy of typically 
50 — 100ms -1 on an absolute velocity scale. The remaining 
minor discrepancies point to how the convection simulations 
can be refined further. 

Key words: Convection - Hydrodynamics - Line: formation - 
Radiative transfer - Sun: granulation - Sun: photosphere 



1. Introduction 

A proper understanding of convection has for a long time been 
one of the greatest challenges for stellar astrophysics. Since 
convection can strongly influence both the stellar evolution and 
the emergent stellar spectra, this uncertainty also carries over to 
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other fields of astrophysics, such as galactic evolution and cos- 
mology. Convection also plays a dominant role in various pro- 
cesses besides energy transport and stellar evolution, perhaps 
most notably mixing of nuclear-processed material, magnetic 
activity, chromospheric heating and stellar pulsations. Convec- 
tion is furthermore of relevance in searches for extra-solar plan- 
ets by introducing differential spectral line shifts that are likely 
to change periodically with solar-like cycles, and may then af- 
fect the measurements of periodic variations induced by orbit- 
ing companions. 

The convection zone in late-type stars extends to the visible 
surface layers and thereby influences the emergent spectrum, 
from which virtually all information about the stars is deduced. 
For the Sun, the directly observable manifestation is granula- 
tion: an evolving pattern of brighter (warmer) rising material 
with darker (cooler) sinking gas in between, driven by the ra- 
diative cooling in a thin surface layer. Convection also leaves 
distinct signatures in the spectral lines in the form of line shifts 
and asymmetries (e.g. Gray 1980; Dravins et al. 1981; Dravins 
1982; Dravins & Nordlund 1990a,b). For late-type stars, the 
bisectors of unblended lines have a characteristic C-shape, 
which results from the differences in line strength, contin- 
uum level, area coverage and Doppler shift between up- (gran- 
ules) and downflows (intergranular lanes). Weak lines normally 
only show the upper part of the C-shape and have convective 
blueshifts (after removing the effects of the gravitational red- 
shift). The cores of stronger lines are formed in or above the 
convective overshoot layers and therefore have smaller or non- 
existent blueshifts (Allende Prieto & Garcia Lopez 1998a,b). 
Solar transition region and coronal emission lines can, how- 
ever, show significant redshifts of several kms -1 (Doschek et 
al. 1976). 

The possibility to obtain spectra of the Sun on an absolute 
wavelength scale (i.e. corrected for the radial velocity of the 
Sun relative to the spectrograph) with exceptionally high spec- 
tral resolution, signal-to-noise ratio (S/N) and wavelength cov- 
erage, makes the Sun an ideal testbench to study the physical 
conditions in the upper parts of the convection zone. For other 
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stars one has previously been restricted to only a few spectral 
lines and poorer quality observations. Fortunately the situation 
has recently improved significantly (e.g. Allende Prieto et al. 
1999, 2000), most notably concerning resolution, S/N and the 
use of carefully wavelength-calibrated spectra covering a large 
portion of the visual region. 

In order to draw the correct conclusions from the observed 
stellar spectrum it is necessary to have a proper understanding 
of how convection modifies the atmospheric structure and the 
spectral line formation. Up to the present, theoretical models 
of stellar atmospheres have been based on several unrealistic 
assumptions, such as restriction to a 1 -dimensional (ID) strat- 
ification in hydrostatic equilibrium with only a rudimentary 
parametrisation of convection through e.g. the mixing length 
theory (Bohm-Vitense 1958) or some close relative thereof 
(Canuto & Mazzitelli 1991). Such models therefore miss en- 
tirely the very nature of convection as an inhomogeneous, non- 
local and time-dependent phenomenon, which is driven by the 
radiative cooling at the surface. Furthermore, spectral synthe- 
sis using classical ID model atmospheres show poor resem- 
blance with observed spectra. In particular, the predicted line 
broadening is much too weak (of course, with no knowledge 
about the convective motions and Doppler shifts, all lines are 
also perfectly symmetric). To partly cure the problem, two ar- 
tificial fudge parameters are introduced which the theoretical 
spectrum is convolved with to mimic the missing broadening: 
micro- and macroturbulence. In this situation, systematic errors 
arising from modeling inadequacies in e.g. chemical analyses 
of stars are far larger than the measuring errors alone. An im- 
proved theoretical foundation for the interpretation of stellar 
spectra is therefore highly desirable, to complement the im- 
pressive recent observational advances in terms of S/N, spectral 
resolution and limiting magnitudes following the new genera- 
tion of large-aperture telescopes, like the ESO VLT and Mc- 
Donald HET, and detectors. 

With the advent of modern supercomputers, self-consistent 
3D radiative hydrodynamical simulations of the surface con- 
vection in stars have now become tractable (e.g. Nordlund & 
Dravins 1990; Stein & Nordlund 1989, 1998; Asplund et al. 
1999; Trampedach et al. 1999). These simulations invoke no 
free parameters like mixing length parameters except the nu- 
merical viscosity for stability purposes, yet succeed in repro- 
ducing key observational diagnostics such as granulation topol- 
ogy and statistics (Stein & Nordlund 1989, 1998), and helio- 
seismic properties (Rosenthal et al. 1999). In the present pa- 
per we present a detailed study of the spectral line formation 
of Fe lines using such solar surface convection simulations as 
3D model atmospheres. The resulting line shapes, shifts and 
asymmetries are found to compare very satisfactory with the 
observations, lending strong support to the realism of the simu- 
lations. The derived solar Fe abundances from weak and strong 
Fe I and Fe II lines are presented in an accompanying article 
(Asplund et al. 2000b, hereafter Paper II). Subsequent papers 
in the series will deal with Si lines (Asplund 2000, hereafter 
Paper III) and observed and predicted spatially resolved solar 
lines (Kiselman & Asplund, in preparation, hereafter Paper IV). 



Additionally, hydrogen line formation will be the topic of a fol- 
lowing article (Asplund et al., in preparation, hereafter Paper 
V). 

2. Surface convection simulations 

The 3D model atmospheres of the solar granulation which 
form the basis for the spectral line calculations presented here 
have been obtained with a 3D, time-dependent, compressible, 
radiative-hydrodynamics code developed to study solar and 
stellar surface convection (e.g. Nordlund & Stein 1990; Stein & 
Nordlund 1989, 1998; Asplund et al. 1999). The hydrodynam- 
ical equations of mass, momentum and energy conservation: 



= —v • Vlnp — V • v 



d\np 
~df 

dv P 1 

— = -v • Vv + g VlnP + -V • a 

at p p 

de P 

ttt = -v • Ve V • v + Q rad + Q visc 

at p 
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coupled to the equation of radiative transfer along a ray (in total 
eight inclined rays): 



d/A/dr A = I\ - S x 



(4) 



are solved on a non-staggered Eulerian mesh with 
200x200x82 gridpoints; simulations with resolutions 
100x100x82, 50x50x82 and 50x50x41 have also been 
computed but have not been utilized here since they are 
hampered by less good agreement between predicted and 
observed line profiles (Asplund et al. 2000a). In these 
equations, p denotes the density, v the velocity, g the grav- 
itational acceleration, P the pressure, e the internal energy, 
a the viscous stress tensor, <5 V i SC the viscous dissipation, 
Qrad = A. Jo K \{I\ " S\)dfld\ the radiative heating/cooling 
rate, I\ the monochromatic intensity, 5a = B\(T) the 
corresponding source function, t\ the optical depth and k\ the 
absorption coefficient. 

The physical dimension of the simulation box corresponds 
to 6.0x6.0x3.8 Mm of which about 1.0 Mm is located above 
continuum optical depth unity. Again, initial solar simulations 
extending only to heigths of 0.6 Mm suffered from problems 
in the predicted line asymmetries for strong lines which was 
traced to the influence of the outer boundary and that non- 
negligible optical depths were present in the line cores al- 
ready at the outermost layers. Though not completely removed 
with the current more extended simulations, the problems have 
largely disappeared as discussed further in Sect. || and ^. The 
depth scale has been optimized to provide the best resolution 
where it is most needed, i.e. in those layers with the steepest 
gradients in terms of dT/dz and d 2 T/dz 2 , which for the Sun 
occurs around the visible surface (which is defined to have a 
geometrical depth z = OMm). Through this procedure, suf- 
ficient depth coverage was also obtained in the important op- 
tically thin layers. The simulation box covers about 13 pres- 
sure scale heights and 1 1 density scale heights in depth, while 
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the horizontal dimension is sufficient to include <; 10 gran- 
ules at any time of the simulation. The spatial derivatives are 
computed using third order splines and the time-integration 
is a third-order leapfrog predictor-corrector scheme (Hyman 
1979; Nordlund & Stein 1990). The code was stabilized using 
a hyper-viscosity diffusion algorithm (Stein & Nordlund 1998) 
with the parameters determined from standard hydrodynamical 
test cases, like the shock tube. In this respect these parameters 
are not freely adjustable parameters for the surface convection 
simulations. It is important to realize that changing the numeri- 
cal resolution in effect also varies the effective viscosity, which 
shows that the convective efficiency and temperature structure 
are independent on the adopted viscosity description (Asplund 
et al. 2000a). Periodic horizontal boundary conditions and open 
transmitting top and bottom boundaries were used in an iden- 
tical fashion to the simulations described by Stein & Nordlund 
(1998). The influence of the boundary conditions on the re- 
sults are negligible. In particular the upper boundary has been 
placed at sufficiently great heights not to disturb the line for- 
mation process except for the cores of very strong lines, and 
the lower boundary is located at large depths to ensure that the 
inflowing gas is isentropic and featureless. No magnetic fields 
or rotation were included in the present calculation. 

In order to obtain a realistic atmospheric structure it is cru- 
cial to correctly describe the internal energy of the gas and the 
energy exchange between radiation and gas. For this purpose 
a state-of-the-art equation-of-state (Mihalas et al. 1988), which 
includes the effects of ionization, excitation and dissociation, 
has been used. Since the solar photosphere is located in the lay- 
ers where the convective energy flux from the interior is trans- 
ferred to radiation, a proper treatment of the 3D radiative trans- 
fer is necessary, which has been included under the assumption 
of LTE and using detailed continuous (Gustafsson et al. 1975 
and subsequent updates) and line (Kurucz 1993) opacities. The 
3D equation of radiative transfer was solved at each timestep 
of the convection simulation for eight inclined rays (2 //-angles 
and 4 (^-angles) using the opacity binning technique (Nordlund 
1982). The four opacity bins correspond to continuum, weak 
lines, intermediate strong lines and strong lines; shorter tests 
with eight instead of four opacity bins revealed only very minor 
differences in the resulting atmospheric structures. The accu- 
racy of the binning procedure was verified throughout the sim- 
ulation at regular intervals by solving the full monochromatic 
radiative transfer (2748 wavelength points) in the 1.5D approx- 
imation (i.e. treating each vertical column separately in the flux 
calculations without considering the influence from neighbor- 
ing columns) and found to always agree within 1 % in emergent 
flux. 

It is noteworthy that the convection simulations contain 
no adjustable free parameters besides those used to character- 
ize the stars: the effective temperature T c s (or equivalently, as 
adopted here, the entropy of the inflowing material at the bot- 
tom boundary), the surface acceleration of gravity \ogg and 
the chemical composition. The entropy at the lower boundary 
was carefully adjusted prior to the simulation started in order to 
reproduce the total solar luminosity. The solar surface gravity 



(log g — 4.437 [cgs]) was employed. The chemical composi- 
tion of the gas were taken from Grevesse & Sauval (1998), in 
particular the Fe abundance used for the equation-of-state and 
opacity calculations was 7.50f} The line opacities were inter- 
polated to the adopted Fe abundance using the standard Kurucz 
ODFs with log ep = 7.67 and non-standard ODFs computed 
with log 6p c = 7.51 and without any He (Kurucz 1997, private 
communication, cf. Trampedach 1997 for details of the proce- 
dure). 

3. Spectral line calculations and atomic data 

From the full solar convection simulation, which covers two 
solar hours, a shorter sequence of 50 min with snapshots stored 
every 30 s was chosen for the subsequent spectral line calcu- 
lations. The adopted snapshots have a temporally averaged ef- 
fective temperature very close to the observed nominal T ff = 
5777 ± 3 K for the Sun (Willson & Hudson 1988): 5767 ± 21 K 
where the quoted range is the standard deviation of the indi- 
vidual snapshots. For our present purposes the time coverage 
is sufficient to obtain statistically significant spatially and tem- 
porally averaged line profiles and shifts, as verified by test cal- 
culations after dividing the snapshots in various sub-groups. 
Even as short sequences as 10 min produced line asymmetries 
and shifts to within 50 m s _1 of the full calculations. In order to 
improve the vertical sampling, the simulation was interpolated 
to a finer depth scale with the same number of depth points 
but only extending down to layers with minimum continuum 
optical depths of r cont = 300 (« 700 km below r con t = 1 
rather than 2.9 Mm in the original simulation). Prior to the line 
transfer computations, the snapshots were also interpolated to 
a coarser grid of 50 x 50 x 82 but with the same horizontal di- 
mension to save computing time; various tests assured that the 
procedure did not introduce any differences in the spatially av- 
eraged line shapes or asymmetries. The Doppler shifts intro- 
duced by the convective motions in the 3D model atmosphere 
were correctly accounted for in the solution of the 3D spec- 
tral line transfer. In most cases, intensity profiles at the center 
of the solar disk (// = 1.0) were considered, which have been 
calculated for every column of the interpolated snapshots, be- 
fore spatial and temporal averaging and normalization. A few 
lines where, however, computed under different viewing angles 
to allow a disk-integration for comparison with published solar 
flux atlases (Kurucz et al. 1984). The assumption of LTE in the 
ionization and excitation balances and for the source function 
in the line transfer calculations (£„ = B u ) have been made 
throughout in the present study. The background continuous 
opacities were calculated using the Uppsala opacity package 
(Gustafsson et al. 1975 and subsequent updates). 

Fe lines are the most natural tools to study line shapes and 
asymmetries in stars: Fe is an abundant element with a complex 
atomic structure which ensures there are many useful lines of 
the appropriate strength; there exists accurate laboratory mea- 

1 On the customary logarithmic abundance scale defined to have 
loge H = 12.00 
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surements for the necessary atomic data such as transition prob- 
abilities and wavelengths (e.g. Blackwell et al. 1995; Nave et al. 
1994); the Fe nuclei have a large atomic mass which minimizes 
the thermal velocities; the influence from isotope and hyper- 
fine splitting should be negligible (e.g. Kurucz 1993); and LTE 
is a reasonable approximation for Fe (e.g. Shchukina & Tru- 
jillo Bueno 2000), at least for ID model atmospheres of solar- 
type stars (3D NLTE studies of Fe may, however, reveal larger 
departures from NLTE as speculated by e.g. Nordlund 1985 
and Kostik et al. 1996). By studying a large sample of neu- 
tral and ionized Fe lines with different atomic data and there- 
fore line strengths, the solar photospheric convection at vary- 
ing atmospheric layers can be probed by analysing the result- 
ing line shifts and asymmetries. The Fe lines and their atomic 
data are the same as described in detail in Paper II. In particu- 
lar, accurate laboratory wavelengths for the Fe I and Fe II lines 
were taken from Nave et al. (1994) and Johansson (1998, pri- 
vate communication), respectively. The final profiles have been 
computed with the individual Fe abundances derived in Paper 
II. 

For comparison with observations, the solar FTS disk- 
center intensity atlas by Brault & Neckel (1987) (see also 
Neckel 1999) has been used, due to its superior quality over the 
older Liege atlas by Delbouille et al. (1973) in terms of wave- 
length calibration (Allende Prieto & Garcia Lopez 1998a,b) 
and continuum tracement. For flux profiles the solar atlas by 
Kurucz et al. (1984) has been used, which is also based on 
FTS-spectra with a similar spectral resolution as the disk-center 
atlas. The wavelengths for the observed profiles have been ad- 
justed to remove the effects of the solar gravitational redshift 
(633 ms -1 ). All spatially averaged theoretical profiles have 
been convolved with an instrumental profile to account for the 
finite spectral resolution of the observed atlas. Since the atlas 
was acquired with a Fourier Transform Spectrograph (FTS), 
the instrumental profile corresponds to a sinc-function with 
A/AA ~ 520 000 in the visual rather than the normal Gaus- 
sian (e.g. Gray 1992). The additional instrumental broadening 
has only a minor (but not entirely negligible) effect on the re- 
sulting line asymmetries due to the high resolving power of the 
FTS. 
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Fig. 1. Spatially resolved profiles and bisectors for the Fe I 
608.2 line. The lines are both stronger and have higher contin- 
uum intensities in the blueshifted granules compared with the 
red-shifted intergranular lanes. The largest vertical velocities 
in the photosphere are encountered in the downflows. The in- 
tensity scale is normalized to the spatially averaged continuum 
level. The thick solid lines correspond to the spatially averaged 
profile and bisector 



4. General features of 3D line formation 



predicted profiles lend strong support to the realism of the sim- 
ulations. 



The convective motions and the atmospheric inhomogeneities 
leave distinct fingerprints in the spectral lines, which can be 
used to decipher the conditions in the line-forming layers pro- 
duced by the convection (e.g. Dravins et al. 1981; Dravins 
& Nordlund 1990a,b). Line strengths of weak lines are pre- 
dominantly determined by the average atmospheric tempera- 
ture structure, while the line widths reflect the amplitude of the 
Doppler shifts introduced by the velocity fields. Line shifts and 
bisectors are created by the correlations between temperature 
and velocity and the details of the convective overshooting, as 
well as the statistical distribution between up- and downflows. 
Therefore, obtaining a good agreement between observed and 



4.1. Spatially resolved profiles 

Although Paper IV and V in the present series of articles 
will discuss in detail observed and predicted spatially resolved 
lines, a brief excursion is still warranted here in order to inter- 
pret the spatially averaged profiles and bisectors presented in 
Sects. |, | and |. 

Spatially resolved profiles take on an astonishing range of 
shapes and shifts spanning several km s _1 , as illustrated in Fig. 
I]. The intensity contrast reversal in the higher layers of the pho- 
tosphere, i.e. granules become dark while intergranular lanes 
become bright a few hundred km above the visible surface, is 
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Fig. 2. The spatially and temporally averaged Fel 608.2 line 
when artificially removing all convective velocities in the simu- 
lations (diamonds). In comparison with the solar intensity atlas 
(solid line, Brault & Neckel 1987) the predicted line is much 
too narrow and lacks the correct line shift and asymmetry. The 
result when including the self-consistent velocity field is shown 
in Fig. | 



Fig. 3. The spatially and temporally averaged Fe I 608.2 bisec- 
tor assuming three different abundances: log ep e = 7.00, 7.50 
and 8.00 (or equivalently with three different loggj-values). 
Note that the upper parts of the bisectors do not coincide for 
the different line strengths since the whole region of line for- 
mation is shifted outwards for stronger lines 



clearly seen in the cores of the resolved profiles. The strengths 
of spatially averaged profiles are normally biased towards the 
granules, since the upflows in general are brighter (high con- 
tinuum intensity), have steeper temperature gradients and have 
a larger area coverage than the downfiows. These trends are 
also observationally confirmed, which suggests that the com- 
bination of 3D hydrodynamical model atmospheres and LTE 
is appropriate for most Fe lines (Kiselman 1998; Kiselman & 
Asplund 2000; Paper IV). 

Spatially resolved bisectors are not at all typical of the spa- 
tially averaged bisectors, which are merely the result of the sta- 
tistical distribution of individual profiles. Rather than the char- 
acteristic C -shape bisectors, individual bisectors normally have 
inverted \-shape bisectors in granules and /-shapes in inter- 
granular lanes, as seen in Fig. [jj This reflects the in general 
increasing vertical velocities with depth in the photosphere. 
However, due to the meandering motion of the downfiows, oc- 
casionally the line-of-sight passes through both upward and 
downward moving material which causes large variations for 
certain columns. 

4.2. Spatially averaged disk-center profiles 

The widths of spatially averaged spectral lines, which clearly 
exceed the natural and thermal broadening, predominantly 
arise from the velocity amplitude of the granules and inter- 
granular lanes and to a lesser extent from photospheric oscil- 
lations. A demonstration of the importance of the non-thermal 
Doppler broadening is presented in Fig. ||, which shows the re- 
sulting spatially averaged profile from the 3D simulations but 
with all convective velocities artificially set equal to zero in 
the line calculation; thereby the predicted profile closely re- 
semble those calculated with classical ID model atmospheres. 
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Fig.4. The bisectors of the Fel 608.2 (weak) and 621.9 nm 
(strong) lines (solid lines) which both have Xoxc = 2.2 eV. 
Also shown are the corresponding (artificial) Fe I lines with 
X = 1.2 eV (dotted lines) and \ = 3.2 eV (dashed lines). The 
fake lines have all other transition properties the same as the 
original lines but with the (//-values (~ ±1.0 dex) adjusted to 
return the same line depths. All lines have here been computed 
for the same 5 min sequence from the full convection simu- 
lation, which means that the bisectors differ slightly from the 
average of the whole 50 min simulation 



Clearly, without the Doppler shifts the line is much too nar- 
row, which requires additional broadening in the form of micro- 
and macroturbulence to be introduced. The poor agreement 
between observations and predictions when not including the 
self-consistent velocity field as shown in Fig. ||, should be con- 
trasted with the excellent fit shown in Fig. |J 
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Fig.5. The bisectors of the Fel 608.2 (weak) and 621. 9nm 
(strong) lines (solid lines). Also shown are the correspond- 
ing (artificial) Fe I lines with A = 500 nm (dotted lines) and 
A = 700 nm (dashed lines). The fake lines have all other tran- 
sition properties the same as the original lines but with the gf- 
values adjusted to return the same line depths. All lines have 
here been computed for the same 5 min sequence as in Fig. Q 

The individual line bisectors depends on the details of the 
line formation and thus on transition properties such as log gf, 
X cxc and A in an intricate way, as illustrated in Fig. ^, Q and |5|. 
When discussing mean bisectors (e.g. Gray 1992; Allende Pri- 
eto et al. 1999; Hamilton & Lester 1999) it is therefore impor- 
tant to consider only appropriate subgroups consisting of lines 
with similar characteristics to avoid introducing errors when 
interpreting the results in terms of convective properties. As a 
corollary, it follows that reconstructing a mean bisector by av- 
eraging bisectors or using shifts of lines of different strengths 
will in general not recover in detail the individual bisectors, 
as exemplified in Fig. [|. As expected, for line depths <; 0.5 
the bisectors closely coincide for lines of different strengths 
due to the disappearance of the influence from the convective 
inhomogeneities. Similarly, decreasing the excitation potential 
tends to shift the line formation outwards, producing less con- 
vective blueshifts for a given line depth (Fig. Q), while decreas- 
ing the wavelength increases the brightness contrast and makes 
the temperature gradients steeper, resulting in more vigorous 
convective line asymmetries (Fig. |]). 

The temporal evolution of the granulation pattern also in- 
troduce changes in the line profiles. Solar granules have typical 
lifetimes of about 10 min but since the numerical box covers 
typically <; 10 granules at any time, the influence of growing 
and decaying granules are relatively modest, as illustrated in 
Fig. |6| Since the emergent T c $ is not enforced but rather is the 
result of the evolving granulation pattern, the continuum level 
varies slightly (« ±2%, corresponding to T c ff. rms = 20 K) 
throughout the simulation, as well as the line strength. The 
bisector shapes are only slightly modified by the granulation, 
though the presence of oscillations in the simulation box shifts 
the bisectors back and forth in a regular fashion with a period 




0.0 Lx_ 

-1 n 



-0.4 -0.2 0.0 
Velocity [krn/s] 



Fig. 6. The temporal evolution of the spatially averaged profile 
and bisector for the Fe I 608.2 line. The 40 profiles shown here 
are taken at one minute intervals from the full solar simulation. 
The intensity scale is normalized to the averaged continuum 
level for the full simulation sequence 



corresponding to the solar 5 min oscillations. These numerical 
oscillations have an amplitude of about ±300 m s _1 in the line- 
forming layers. The oscillations additionally broaden the line, 
though without altering the line strength since essentially the 
whole photosphere oscillates in unison without modifying the 
atmospheric structure. 

4.3. Spatially averaged off-center intensity and disk-integrated 
flux profiles 

Although not the main emphasis of the present paper a few 
Fe lines have been computed at different viewing angles (4 
/j-angles and 4 (^-angles) in order to enable a disk-integration 
to obtain flux profiles. A disk-integration requires further that 
the rotational broadening (1.8 km s -1 in the case of the Sun) is 
taken into account. The line formation of flux profiles is more 
complex than for intensity profiles due to the contributions 
from different disk positions. The well-known limb-darkening 



M. Asplund et al.: Solar line formation: I. Fe line shapes, shifts and asymmetries 



7 



1.0 

>- 0.8 











0.96 






f. = 


0.82 ~ 








0.57 




F=l \ 








A = 608.271 [nm]\ 


h = 






X = 2.223 [eV] \. * 








log gf = -3.573 




4 




log E = 7.50 




1 





608.250 608.260 608.270 608.280 608.290 
Wavelength [nm] 

1.2 



1 .0 - 




0.2 - 

o.o r " 

-0.50 -0.40 -0.30 -0.20 -0.10 0.00 0.10 0.20 
Velocity [km/s] 

Fig. 7. The spatially and temporally averaged Fel 608.2 pro- 
file and bisector at different viewing angles. The line has been 
computed at 4 /i-angles and 4 (^-angles. The profiles have been 
normalized to the continuum intensity at disk center. The thick 
solid lines represent the profile and bisector of the line at disk- 
center. Note the "limb effect" in the bisectors: smaller convec- 
tive blue-shift toward the limb 

decreases the continuum intensities towards the limb but the 
lines also tend to be weaker due to the more shallow tempera- 
ture gradient in the upper atmosphere. Furthermore, the granu- 
lation contrast decreases towards the limb, since higher layers 
with progressively smaller influence of the granulation is seen. 
At a given instant the line asymmetries vary significantly more 
towards the limb as the line-of-sight may pass through both 
granules and intergranular lanes and also the horizontal veloc- 
ities may introduce Doppler shifts. Thus the Doppler shifts are 
less well correlated with the background continuum intensity, 
which introduces a more random nature of the resulting bisec- 
tors for inclined line-of-sights. 

Fig. ^ illustrates the variation of the spatially and tempo- 
rally averaged intensity profiles and bisectors of a weak Fe I 
line at various viewing angles, which make up the necessary 
profiles for a disk-integration. Also seen is the limb-effect 
(Halm 1907): the wavelength of solar lines increases towards 



the limb such that at small /i the typical convective blueshift 
seen at disk-center has disappeared completely or even been re- 
versed into a small red-shift when removing the effects of grav- 
itational redshift and rotation. The limb-effect is the same ef- 
fect as present between strong and weak lines, namely that the 
core of the lines are formed in high enough photospheric lay- 
ers where the granulation contrast and velocities have largely 
vanished. The red-shifted cores at very small fi, which are 
both observed and predicted with the 3D model atmospheres, 
are likely due to a bias for receding horizontal velocities to 
be viewed against the higher temperatures above intergranu- 
lar lanes while the approaching gas will preferentially be seen 
against the lower temperatures above granules due to the tem- 
perature reversal in the higher photospheric layers (Balthasar 
1985). 



5. Solar Fe line shapes 

The Doppler shifts introduced by the convective flow veloc- 
ities in the photosphere cause significant line broadening be- 
yond the thermal, radiative and collisional broadening. Fig. || 
shows a few examples of spatially and temporally averaged 
Fe I and Fe II lines at disk-center calculated using the 3D solar 
simulation together with the corresponding observed intensity 
profiles; additional examples are given in Paper II. Clearly the 
agreement is very satisfactory for unblended Fe I and Fe II lines 
with almost a perfect match. In contrast, theoretical ID profiles 
are clearly discrepant in spite of the presence of both a micro- 
and macroturbulence (as exemplified by the Fel 621.9 nm line 
which is shown in the middle panel of Fig. ||). The residual in- 
tensities amount to only < 1%, which is accomplished without 
any micro- or macroturbulent broadening. The FWHM of the 
Fe lines typically agree to within 1%. Fig. ^ shows the disk- 
integrated flux profiles of two Fe I lines as compared with the 
solar flux atlas of Kurucz et al. (1984) and again the agreement 
is very encouraging. In fact, given the good agreement in gen- 
eral for both intensity and flux profiles, blending lines, which 
may otherwise have gone undetected in a ID analysis, are eas- 
ily identified. 

Although the overall agreement is very satisfactory, it is 
clear from a closer inspection of Figs. ^| and^ that there are sys- 
tematic discrepancies in the line profiles which are appearent in 
most intensity and flux profiles, in particular the weaker lines. 
The cores of the predicted lines tend to be slightly too shallow 
while the near line wings are somewhat too broad, which sug- 
gests a slightly over-estimated rms vertical velocity amplitude 
in the solar simulation. The slightly problematic line cores may 
also signal departures from LTE (cf. Rutten & Kostik 1982), 
which is more likely to affect the cores than the wings. Fur- 
thermore, the cores of intermediate strong lines tend to be dis- 
placed compared with observations; the latter feature will be 
discussed further in Sects. |6| and [7| It should not come as a 
surprise that the cores of the stronger lines show discrepancies, 
since the highest atmospheric layers are likely the least realistic 
due to still missing ingredients in the simulations and spectral 
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Fig. 9. The predicted flux profiles of the Fel 608.2 and 
621.9 nm lines (diamonds) compared with the solar flux atlas 
(solid lines, Kurucz et al. 1984). The theoretical lines have been 
disk-integrated using a solar rotational velocity of ii ro tsim = 
1.8 km s -1 and convolved with sinc-function to account for the 
finite spectral resolution of the FTS-atlas 

synthesis in terms of e.g. departures from LTE and the inexact 
line blanketing treatment in the actual convection simulation. 

The minor disagreements shown by the Fe lines are very 
important, since they point to how the simulations can be 
improved further. Prior to the convection simulations used 
here, we have carried out several similar solar simulations 
which differed from the present ones, most notably in terms 
of equation-of-state and opacities (Gustafsson et al. 1975 with 
subsequent updates vs. Mihalas et al. 1988 and Kurucz 1993), 
height extension (z m ; n = —0.6 vs. —1.0 Mm) and numer- 
ical resolution (100x100x82, 50x50x82 and 50x50x63 
vs. 200x200x82). In all cases they suffered from more pro- 
nounced problems in terms of line asymmetries, which were 
subsequently addressed with the more refined and improved 
simulations. However, it is noteworthy that the overall shapes 
of weak Fe lines illustrated in Figs. || and ^ were even better 
described with the previous equation-of-state and line opacities 
(Gustafsson et al. 1975 with subsequent updates), as seen Fig. 
[io| , although with a slightly smaller Fe abundance. One may 
speculate that the older equation-of-state better describes the 
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Fig. 10. The predicted spatially and temporally averaged profile 
of the Fe 1 608.2 line (diamonds) using a solar simulation based 
on a previous equation-of-state and line opacities (Gustafsson 
et al. 1975 with subsequent updates instead of Mihalas et al. 
1988 and Kurucz 1993) compared with the solar disk-center 
intensity atlas (solid lines, Brault & Neckel 1987). The agree- 
ment is even better than the corresponding profile shown in Fig. 



conditions typical of the line-forming layers than the more re- 
cent version (Mihalas et al. 1988) which has been optimized for 
stellar interiors. But it may also be coincidental, since the dif- 
ferences in line profiles are relatively small and there are also 
additional differences in terms of resolution (253 x 253 x 163) 
and z m i n (—0.6 Mm) between the two simulations. A further 
investigation into the matter would, however, be interesting. 



6. Solar Fe line shifts 

A major advantage with solar observations compared with cor- 
responding stellar observations is the existence of high-quality 
spectral atlases given on an absolute velocity scale, which is 
possible since the differential radial velocity between the Sun 
and the Earth can be accurately corrected for (e.g. Kurucz et al. 
1984; Brault & Neckel 1987; Neckel 1999). Furthermore, the 
well-determined solar mass and radius allow the solar grav- 
itational redshift of 636 ms -1 , or 633 ms -1 for light inter- 
cepted on Earth (Lindegren et al. 1999), to be estimated, leav- 
ing the remaining shifts to be attributed to convection as pres- 
sure shifts and similar line shifts are of much lesser importance 
(e.g. Allende Prieto 1998). Solar spectral lines show different 
line shifts depending on the typical depths of formation of the 
line cores. Unfortunately, the remaining uncertainties in the in- 
dividual laboratory wavelengths and possible blends cause a 
scatter of about 100 m s _1 in the observed line shifts, in partic- 
ular for the weaker lines. 

Through the self-consistently calculated convective flows 
in the solar simulations the predicted line shifts can be di- 
rectly compared with observations on an absolute wavelength 
scale. Figs. [H] and |l2| show the calculated line shifts for Fe I 
and Fe II lines, respectively; we prefer to plot the shifts vs 



M. Asplund et al.: Solar line formation: I. Fe line shapes, shifts and asymmetries 



9 



0.0 



-0.4 



.= -0.6 



Iff o*<js n o$ m 



Equivalent width [pm] 



1 



0.0 



-0.4 - 



S -0.6 




Equivalent width [pm] 



E 0.2 



-0.4 



* . 4 



Equivalent width [pnn] 



Fig. 11. Upper panel: The predicted (open circles) and ob- 
served (stars) line shifts for Fe I lines as a function of equiv- 
alent widths. The equivalent widths are taken from Blackwell 
et al. (1995), Holweger et al. (1995) and Moore et al. (1966), in 
this order of preference. Lower panel: The differences between 
predicted and observed line shifts. The gravitational red-shift 
of 633 m s -1 has been subtracted from the observed line shifts. 
The agreement is very good for weak and intermediate strong 
lines but becomes progressively worse for stronger lines. Most 
of the scatter for the weaker lines is probably due to uncertain- 
ties in laboratory wavelengths and blends 



line strengths rather than vs line depths as Hamilton & Lester 
(1999), since due to saturation the latter method will not de- 
tect clearly the gravitional redshift plateau (Allende Prieto & 
Garcia Lopez 1998a). Weaker lines show as expected more 
prominent blueshifts and Fe II lines more so than Fe I lines due 
to their in general deeper layers of formation. The maximum 
predicted blueshift for Fe I lines is 550 ms -1 while it reaches 
800 ms -1 for Fe II lines. A trend between line shift and excita- 
tion potential is present (Fig. |l3|) but much less pronounced due 
to the obscuration introduced by the W\- and A-dependencies 
(cf. Fig. [| | and Hamilton & Lester 1999). 

As clear from Figs. [H] and |l2]the predicted line shifts agree 
well with observations for weak and intermediate strong lines. 
The average difference for weak (W\ Si 6 pm) Fel lines is 
0±53m s -1 (0.000±0.053kms -1 ), while when including also 
the intermediate strong lines (W\ 10pm) the corresponding 
difference increases to 51 ± 81ms -1 ; in these estimates we 
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Fig. 12. Same as Fig. |TT] but for the 15 Fe II lines. The equiva- 
lent widths are taken from Hannaford et al. (1992) 



have not included the Fe 1 666.8 nm line (W\ = 0.89 pm) as we 
suspect that it has an erroneous laboratory wavelength, since 
it is significantly more discrepant (300 ms -1 ) than other lines 
with similar strengths. The situation is slightly worse for the 
Fe II lines with a mean difference of— 64±85ms -1 for all 15 
lines with W\ ^ 10 pm. However, the laboratory wavelengths 
for Fe II lines are likely of somewhat poorer quality than the 
Fe II lines and there are unfortunately only few available lines 
of the appropriate strengths. 

As evident from Figs, [n] and [j~2|, the correspondance 
between predictions and observations becomes progressively 
worse for stronger Fe I and Fe II lines; a trend with W\ in the 
line shift differences may also extend to weaker Fe II lines, ac- 
cording to Fig. 12 We attribute this effect predominantly to the 



influence of the outer boundary in the spectral line calculations 
but departures from LTE may also play a role. The cores of 
these strong lines have significant optical depths already at the 
uppermost depth layers where the average vertical velocity is 
directed downwards, as shown in Fig. [l4]. This follows from 
having on average no net mass flux, since the downward mov- 
ing material in general has smaller densities and thus larger ve- 
locities in these layers. Since the temperature-velocity correla- 
tions present further in has essentially completely disappeared 
at these large heights, this results in convective redshifts, which 
do not seem to be present in the observational evidence (Al- 
lende Prieto & Garcia Lopez 1998a). The discrepancy was not- 
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Fig. 13. Same as Fig. |ll]but as a function of excitation potential 
of the lower level. In general high excitation lines are formed 
in deeper layers and thus have larger convective blueshifts. The 
trend is not as striking as in Fig. [Fy, since also e.g. the oscillator 
strength determines the line-formation depth 



icably larger in the simulations with smaller extension (maxi- 
mum height of 0.6 Mm) and numerical resolution (Asplund et 
al. 2000a). We therefore consider the effect as a numerical ar- 
tifact, which should be possible to further reduce by using an 
improved treatment of the outer boundary and higher numerical 
resolution. In spite of the remaining shortcomings, differential 
line shifts are very accurate; the main uncertainty is of obser- 
vational (blends and laboratory rest wavelengths) nature. 

7. Solar Fe line asymmetries 

Line asymmetries carry additional information to line shifts, 
since the lines trace the atmospheric structure throughout the 
line forming region and not only the higher layers where the 
core is formed and thus the line shift. Fig. [D] shows a few ex- 
amples of the predicted line asymmetries and the correspond- 
ing observed bisectors. It should be noted that both types of 
bisectors are on an absolute wavelength scale and the synthetic 
bisectors have therefore not been shifted in velocity to match 
the observations. In order to achieve such a remarkable agree- 
ment it is necessary to have both a very accurate description of 
the atmospheric structure and the details of line formation as 
well as very high quality laboratory wavelengths. Clearly the 



Fig. 14. The time averaged horizontal average (solid) and rms 
(dashed) vertical velocity in the solar simulations. Positive ver- 
tical velocities indicate downward motion. Note the upturn in 
(v z ) which reaches 400 m s ~ 1 in the uppermost parts of the at- 
mosphere. The horizontal part of the curves at the upper bound- 
ary is a consequence of the boundary formulation, which spec- 
ifies the same vertical velocity in the two top layers 



result is very satisfactory. The excellent correspondance in Fig. 
[l5| is not only fortuitous, as is apparent from an inspection of 
Fig. which shows the differences in observed and predicted 
line asymmetries for all the 67 computed Fe I lines. Under ideal 
conditions the differences should all be vertical lines with no 
velocity offset, which in fact is not far from the truth. In partic- 
ular the weaker lines show excellent agreement, while the situ- 
ation becomes progressively worse for the cores of the stronger 
lines, reflecting the shortcomings already discussed in Sect. |[ 
This discrepancy also affects the bisectors closer to the contin- 
uum for the stronger lines, causing the bisector differences to 
be predominantly positive. The situation for the Fe II lines are 
shown in Fig. |l7[ which again is very satisfactory. Both with an 
inferior resolution and height extension in the convection sim- 
ulations the resulting bisectors are of noticably lower quality 
when comparing with observations (Asplund et al. 2000a). The 
good overall agreement in terms of line asymmetries therefore 
lend very strong support to the realism of the convection simu- 
lations. 

In fact it is very easy to detect problematic lines due to 
erroneous laboratory wavelengths (large velocity offset) and 
blends (discrepant bisector shape) when comparing line asym- 
metries. It can be mentioned that initially the wavelengths 
for the Fel 697.2 and 718.9 nm lines were accidentally taken 
from the VALD database rather than from Nave et al. (1994), 
which differed by only 16 and 5 mA (corresponding to 0.7 and 
0.2kms _1 ), respectively, though immediately detected when 
analysing the line asymmetries. Furthermore, the Fe II wave- 
lengths provided by Johansson (1998, private communication) 
were found to be of significantly higher quality than the ones 
given in Hannaford et al. (1992). 
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Fig. 15. The predicted (solid line) and observed (line with error 
bars) bisectors of the Fe I 680.4, 627.1 and 624.0 nm (in order 
of increasing line strength) lines. It should be emphasized that 
the velocity scale is absolute for both the computed and mea- 
sured line asymmetries, i.e. no arbitrary wavelength shifts have 
been applied in order to bring the two into agreement. Clearly 
the predicted bisectors agree rather satisfactory with the obser- 
vations. In fact, the close resemblance for these lines is partly 
fortuitous since both the laboratory wavelengths and the wave- 
length calibration of the solar intensity atlas may only be accu- 
rate to about 30 — 50 ms -1 . In comparison, classical ID model 
atmospheres will of course only produce vertical bisectors at 
zero absolute velocity 

8. The nature of macro- and microturbulence 

The concepts of micro- and macroturbulence are introduced in 
ID analyses in order to account for the missing line broadening 
on length scales less than and larger than a unit optical depth, 
respectively. Naturally such a simplistic division is artificial 
since the motions occur on a range of scales. Furthermore, they 
are supposed to represent turbulent motions and thus are nor- 
mally assumed to be isotropic. In reality the appearance of the 
photospheric granular velocities is clearly more laminar than 
turbulent with distinct upflows and downflows, a direct con- 
sequence of the strong density stratification (Nordlund et al. 
1997). Given the evidences presented here and in Papers II and 
III, there appears to be no need to invoke any macro- or micro- 
turbulence in spectral syntheses based on realistic 3D model 
atmospheres. 

The excellent agreement between predicted and observed 
line profiles shown in Fig. || implies that the theoretical lines 
have the correct widths without the use of any macroturbu- 
lence. The classical concept of macroturbulence can therefore 
be fully explained by the self-consistently calculated convec- 
tive velocity fields and the stellar oscillations. It is less obvious 
that the simulations properly account for the small-scale mo- 
tions normally referred to as microturbulence given the finite 
numerical resolution which may or may not resolve all signif- 
icant velocities. We believe, however, that the currently best 



Fig. 16. Differences between predicted and observed bisectors 
for the Fe I lines used for the present analysis; a gravitational 
red-shift of 633 m s _1 has been subtracted for the observed bi- 
sectors. Except for a tendency for the predicted strong lines 
(line depths greater than 0.5) to have slightly too little blueshift 
or too much redshift, as also obvious in Fig. [TIL the agreement 
is clearly very satisfactory. The line furthest to the left is Fe I 
666.8 nm, whose laboratory wavelength therefore can be sus- 
pected to be slightly in error 
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Fig. 17. Same as Fig. [l^ but for the 15 Fe II lines 

solar simulations described here are of sufficiently high resolu- 
tion to describe also the most important effects of these small- 
scale motions. Firstly, solar simulations with different numeri- 
cal grid resolutions indicate that the velocity distributions have 
essentially converged already at a resolution of 200 x 200 x 82, 
which is also apparent from the predicted line shapes (Asplund 
et al. 2000a). The contribution from unresolved scales should 
therefore be very minor. We emphasize that a high resolution is 
needed in the construction of the 3D model atmospheres to al- 
low the high-velocity tails of the velocity distributions, but that 
smaller resolutions can be used for the spectral synthesis, as 
also verified by extensive testing of the line formation in model 
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atmospheres interpolated to various resolutions from the orig- 
inal 200 x 200 x 82 data cubes prior to the analysis presented 
here. Secondly, the derived Fe II and Si abundances show no 
trend with line strength when using profile fitting (Papers II and 
III), which suggests that the velocities important for the line 
broadening have already been accounted for without resorting 
to the use of extra microturbulent-like velocities. There is, how- 
ever, a minor trend for Fe I lines (Paper II) but since Fe II and 
Si lines should be affected similarly yet show no dependence 
with line strength, we attribute the problem with the Fe I lines 
rather to signatures of departures from LTE. Fe I lines should 
be more susceptible to NLTE effects while Fe II lines are essen- 
tially immune to such effects for solar-type stars (Shchukina & 
Trujillo Bueno 2000). From the results presented here and in 
Papers II and III there is no indication that any extra micro- 
turbulent broadening must be included in a fashion similar to 
that of Atroshchenko & Gadun (1994). We attribute this differ- 
ence mainly to our use of much higher resolution (their solar 
simulations had only 30 3 or 32 3 grid-points), height extension 
(their spectral line calculations only extend up to about 400 km 
above tr OS s — 1) and temporal coverage (their calculations 
were restricted to only one or two snapshots), which allow our 
simulations to better describe the full effects of the convective 
motions. Therefore, either from a pragmatic point of view or 
by advocating the principle of Occam's razor, both macro- and 
microturbulence appear redundant in 3D analyses. 

Thus the predominant explanation is the same for micro- 
and macroturbulence, namely the photospheric granular veloc- 
ity field and temperature inhomogeneities. Additionally, pho- 
tospheric oscillatory motions play a role in the overall macro- 
turbulent broadening of the lines without affecting the line 
strengths. The main component for the microturbulence there- 
fore does not at all arise from the microscopic turbulent mo- 
tions but rather from (gradients in) convective motions that 
are resolved with the current simulations. The small scale en- 
ergy cascades and turbulence associated with a high Reynolds- 
number plasma are present in the solar convection zone but 
their intensities are very small in the upflows to which the line 



strengths are strongly biased (Sect. 4. 1 ) and therefore they do 



not influence the line formation significantly. 
9. Concluding remarks 

The good agreement with observed line shapes, shifts and 
asymmetries, lend very strong support to the realism of the 
3D convection simulations. No doubt the 3D predictions are 
superior to those obtained from ID analyses. In particular the 
use of mixing length parameters, equivalent widths, macro- and 
microturbulence no longer appear to be needed. Therefore de- 
rived results, such as elemental abundances, should be more re- 
liable. In spite of the minor remaining shortcomings, the overall 
significant accomplishment is therefore still obvious: starting 
from only the well-known radiative-hydrodynamical conserva- 
tion equations and with no adjustable free parameters besides 
the treatment of the numerical viscosity in the construction of 
the 3D model atmospheres, detailed line profiles and asymme- 



tries can be predicted which agree almost perfectly with obser- 
vations and furthermore are far superior to classical ID pre- 
dictions with several tunable parameters. It should be stressed, 
however, that the numerical viscosity is merely introduced for 
numerical stability purposes and is determined from standard 
hydrodynamical test cases with no adjustments allowed to im- 
prove the agreement with observations. In this respect the vis- 
cosity is not a freely adjustable parameter like e.g. the various 
mixing length parameters in ID models. 

It is important to emphasize that this accomplishment is 
only possible if the convection simulations are highly realistic, 
both in terms of input physics (equation-of-state, opacities etc) 
and numerical details (numerical and physical resolution, ex- 
tension, boundary conditions, radiative transfer treatment etc). 
From our various experiments and test calculations we can con- 
clude that of special importance is the dimension (2D is not ad- 
equate for spectral synthesis, Asplund et al. 2000a), resolution 
(even ~ 100 3 is not quite sufficient for detailed line shapes, As- 
plund et al. 2000a) and height extension (to limit the influence 
from the outer boundary) of the numerical box. Furthermore, 
in order to achieve the correct temperature structure it is im- 
portant to include the effects of line -blanketing in the 3D radia- 
tive transfer during the simulations (Stein & Nordlund 1998). 
We believe that we have now addressed all of these specific is- 
sues with the new generation of convection simulations, which 
is supported by the very close resemblance with observed line 
profiles, shifts and asymmetries, as presented here. 

All of the above-mentioned features are currently afford- 
able with present-day supercomputers. The obvious disadvan- 
tage with this 3D procedure is of course the time-consuming 
task to perform the necessary 3D convection simulations and 
spectral line calculations, even with the simplifying assump- 
tions of opacity binning and LTE. Furthermore, still only a rel- 
atively small part of the Hertzsprung-Russell diagram of solar- 
type stars has been explored. The situation will, however, im- 
prove significantly during the coming years due to faster com- 
puters and more efficient numerical algorithms. Of particular 
interest will be to extend the modeling to additional metal-poor 
stars (cf. Asplund et al. 1999), A-F type stars and red giants, 
where the granulation is expected to be much more vigorous 
than for the Sun, which should therefore influence the line for- 
mation more. 
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Fig. 8. Some examples of spatially and temporally averaged disk-center Fe I and Fe II lines (diamonds) compared with the ob- 
served solar atlas (solid lines, Brault & Neckel 1987). Also shown are the residual intensities (observed - predicted) to emphasize 
the remaining minor differences; minor blends not included in the spectral synthesis are clearly seen in some of the line wings. 
The predicted profiles have been convolved with a sinc-function to account for the finite spectral resolution of the solar atlas. 
Minor corrections to the observed wavelengths (< 100 ms -1 ) and continuum level (< 1%) have been allowed to improve the 
fits. The Fe lines are almost perfectly matched by the theoretical profiles, which implies that the rms velocity amplitudes in the 
simulations are very close to the real values in the solar photosphere. For comparison the best fit ID profile for the Fe I 621.9 nm 
line is also shown (green solid line in the middle right panel), which has been computed with the Holweger-Miiller (1974) model 
atmosphere, log epd = 7.60 (to achieve the same equivalent width as the 3D profile), £turb = 0.845km s _1 and convolved with 
a Gaussian macroturbulence of 1.6 km s -1 to have the correct line depth (radial-tangential macroturbulence is not applicable for 
intensity profiles, Gray 1992); clearly the agreement is much inferior in spite of the adjustable broadening parameters due to the 
neglect of convective velocities and spatial inhomogeneities 



